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When p- dimensional branes annihilate with antibranes in the early universe, as in brane-antibrane 
inflation, stable (p— 2)-dimensional branes can appear in the final state. We reexamine the possibility 
that one of these (p — 2)-branes could be our universe. In the low energy effective theory, the final 
state branes are cosmic string defects of the complex tachyon field which describes the instability of 
the initial state. We quantify the dynamics of formation of these vortices. This information is then 
used to estimate the production of massless gauge bosons on the final branes, due to their coupling 
to the time-dependent tachyon background, which would provide a mechanism for reheating after 
inflation. We improve upon previous estimates indicating that this can be an efficient reheating 
mechanism for observers on the brane. 

PACS numbers: 98.80.Cq 
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I. INTRODUCTION 



In the last few years significant progress has been 
made in constructing string theoretic cosmological mod- 
els where inflation is driven by the naturally occuriii; 
potentials between D-branes and their antibranes p| 
The formation of lower-dimensional branes at the end of 
inflation can lead to interesting signals: cosmic-string- 
like or higher dimensional defects could be observable 
remnants |4J , possibly providing a rare clue to the stringy 
origin of inflation. A more radical idea was explored in 
H : perhaps our own observable universe is such a defect 
in the higher-dimensional spacetime predicted by string 
theory. Since the stable branes in Type IIB string theory 
have spatial dimensionalities which are odd, a 3-brane 
would have descended from annihilation of 5-branes in 
this picture, and thus they would be codimension-two 
defects in of the effective 6D theory. Codimension-two 
braneworlds have attracted interest lately because of 
their novel features, which might have some bearing on 
the cosmological constant problem 0, Q ■ 

Our interest in this scenario is motivated by questions 
about the efficiency of reheating in brane-antibrane in- 
flation 8]. It is possible that the energy density liber- 
ated from the brane collisions will be converted mostly 
into closed string states, ultimately gravitons, and not 
necessarily into visible radiation [jj]. A generic mecha- 
nism which could avoid this problem was proposed in 
P|, wherein the reheating in D-brane driven inflation is 
due to the coupling of massless gauge fields to a time- 
dependent tachyon condensate, which describes the an- 
nihilation process. However, ref. considered only the 
formation of a tachyon kink instead of the more realistic 
case of a vortex, and it used a somewhat crude ansatz 
for the background tachyon field. The problem of finding 
the actual tachyon background predicted by string theory 
was studied numerically in 01 but no attempt was made 
to improve on the reheating computation. In this paper 



we aim to analytically determine the dynamics of for- 
mation of lower dimensional branes described as tachyon 
defects — both kinks and vortices — and to improve on the 
reheating calculation of p| . 

Let us begin by describing the scenario we have in 
mind. In the simplest version of D-brane inflation a par- 
allel brane and antibrane begin with some separation be- 
tween them in one of the extra dimensions. Although 
parallel branes are supersymmetric and have no force 
between them, the brane-antibrane system breaks super- 
symmetry so that there is an attractive force and hence 
a nonvanishing potential energy. It is the latter which 
drives inflation. Once the branes have reached a criti- 
cal separation one of the stretched string modes between 
the branes, T, becomes tachyonic and the branes become 
unstable to annihilation. The tachyon field starts from 
the unstable maximum T = and rolls towards the vac- 
uum T — > ±oo. However, topological defects may form 
through the Kibble mechanism |j] so that T = stays 
fixed at the core of the defect. These defects are known to 
be consistent descriptions of branes whose dimension is 
lower than that of the original branes For exam- 

ple, the brane-antibrane system has a complex tachyon 
field, leading to vortices which represent codimension-two 
branes. On the other hand, an unstable brane has a real 
tachyon which leads to kinks representing codimension- 
one branes. 

The formation of tachyon defects at the endpoint of D- 
brane inflation is a dynamical process where the tachyon 
couples to gauge fields which will be localised on the de- 
scendant brane. It is thus expected that some radiation 
will be produced by the rolling of the tachyon and the 
problem of reheating becomes quantitative: can this ef- 
fect be efficient enough to strongly deplete the energy 
density of the tachyon fluid so the the universe starts out 
being dominated by radiation rather than cold dark mat- 
ter? It is important to stress that though the situation 
is somewhat analogous to that of hybrid inflation (where 
the tachyon plays the role of the unstable direction in 
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field space which allows for inflation to end quickly) the 
mechanism for reheating is qualitatively different. The 
difference is that in the low energy effective field theory 
which describes the tachyon T, the potential is minimized 
at T = ±00 and there are no oscillations about the mini- 
mum of the potential. In a normal hybrid inflation model, 
T would have a minimum at some finite value and the 
oscillations of T around its minimum would give rise to 
reheating in the usual way. In the present case, the time 
dependence of the background is monotonic, not oscil- 
latory. Reheating thus might seem to resemble gravita- 
tional particle production |l3| rather than the standard 
picture in which the inflaton decays. However, in this 
work we highlight an important difference between re- 
heating through tachyon condensation and gravitational 
particle production, which can make the former much 
more efficient: there is a divergence in the stress-energy 
tensor of the tachyon field within a finite time, which 
corresponds to the formation of the lower-dimension D- 
brane. 

In this paper we study analytically the dynamical for- 
mation of the tachyon vortex and improve the reheating 
calculation, using a slightly simplified model of particle 
production by tachyon condensation, which captures the 
essential physics revealed by the analysis of vortex for- 
mation. In section II we review the formation of tachyon 
kinks which describes the condensation of a brane to a 
brane of codimension-one. In sections III-V we study 
analytically the formation of a tachyon vortex on the 
brane-antibrane pair. Section III introduces Sen's ac- 
tion for the complex tachyon field describing this situ- 
ation. Section IV presents new analytic results for the 
time-dependent, complex tachyon field representing vor- 
tex formation, both near to and far from the vortex core. 
In section V we show that the stress-energy tensor for 
the system splits into a localized, singular piece describ- 
ing the descendant branes, plus a bulk contribution that 
describes the rolling tachyon condensate. Section VI in- 
troduces the effective action for U(l) gauge bosons which 
become localized on the final-state 3-brane, in the rolling 
tachyon background. This provides a model for the vis- 
ible radiation produced during reheating. In section VII 
we calculate the energy density of this produced radia- 
tion on the 3-brane, using some reasonable simplifying 
assumptions. Section VIII gives our conclusions, includ- 
ing speculation about how the final brane-antibrane sys- 
tem could be stabilized. 



II. DYNAMICAL TACHYON KINK 
FORMATION ON UNSTABLE DP-BRANES 



vortex since many of the results are quite analogous. 



A. Effective Field Theory and Equations of Motion 

We will work with the effective action for the tachyon 
on an unstable Dp-brane 0, 1 

S = ~J V{T)y/- det \timn + d M Td N T\ <F +1 x (1) 

where we have set the gauge fields and transverse scalars 
to zero. We use the potential V(T) — t p exp (— T 2 /a 2 ) 
where r p is the tension of a Dp-brane and a — 1\rnol '. 
The value of the constant a is chosen so that the potential 
satisfies the normalization condition 



+00 



V(y) dy = = r p _i 



(2) 



proposed in ^lj • This normalization was used in |Tl| to 
fix the tension of the singular static kink solution of the 
action to correspond to the tension of a D(p — 1)- 
brane. For a time-dependent kink solution we take T to 
be a function of x M = (t, x) so that the action (JIJ becomes 



S = - J V{T)yJl + d^Td^Td v+1 3 



(3) 



Static solutions of the theory ® are well studied in the 
literature 0, 0|. Inhomogeneous solutions have also 
been studied in some detail p| 0, 0, 0. The energy 
momentum tensor for (0) is 

V(T) 



V 1 + dpTdPT 
and the Euler-Lagrange equation of motion is 

V'(T) 



d d Td^Td^T 
dfduT - ' 1 v 



1 



dPTd p T 



V(T) 



= 



(5) 



where V'{T) = ^p-- It is worth noting, as in [T^, that 
the equation of motion (JSJ) is equivalent to conservation 
of energy dpT^ — for nonconstant T since dpT^ v — 

It will be useful in the ensuing 
V(T) 



Qurp 



d, 



analysis to define 



y/l + d^TdPT 



(6) 



In this section we review the dynamical formation of 
a D(p — l)-brane through tachyon condensation on an 
unstable Dp-brane, and derive a few new results. The 
equations of motion in this case are simpler than in the 
case of the vortex and we will use the analysis of this 
section to reinforce our conclusions when we analyze the 



The convention for indices is that upper case roman indices 
{M, N} run over the full space-time coordinates {0, 1, •••,£>}, 
greek indices {fi, u} run over the defect coordinates {0, 1} and 
"hatted" greek indices {(1, 0} run over the remaining spatial coor- 
dinates {2, 3, ■ ■ ■ ,p}. We use metric signature diag(— 1, 1, 1, • ■ ■). 



3 



B. Solutions Near the Core of the Defect 

At the core of the kink we expect the field to stay 
pinned at T = 0. Consider initial data T(t = 0, x) = 
Ti{x) and f(t = 0,x) = f^x) = 2 . One expects the 
field to start to roll where Tj(x) ^ due to the small 
displacement from the unstable maximum V'(T) = 0. 
At t = the equation of motion (JjjJ is 

fi(x)(l + T>{xf) = T!'(x) + 2a^(x)(l + l?(x) 2 ). 

Clearly any point Xq where T^xq) = T/'(x ) = will be 
a fixed point where T(t,xo) = = T(t,Xo) throughout 
the evolution. We restrict ourselves only to considering 
intial data such that sgn(Tj(x)) = sgn(Tj(x)) for all x to 
ensure that the solutions are increasing. 

At the site of the kink (which we take to be xq = 
0) we have T — 0; hence there should always be some 
neighbourhood of the point x = where we can take 
V'(T) ^ so that © yields 



f (1 + T' 2 ) = (1 - T 2 )T" + 2TT'T'. 
This has an increasing solution with T" = 



T(x, t) = x tan — (t — t c ) + — 



(7) 



(8) 



Near the site of the kink the slope of the tachyon field 
diverges as (t c — t)~ x as t approaches the critical time t c , 
similar to the solutions in [lf| . 

The finite-time slope diver genc e was observed both nu- 
merically and analytically in |l0( and leads to the forma- 
tion of a singularity in the energy density at t — t c . This 
effect was also found in an exact string theoretic calcu- 
lation in J2jj. As t — > t c we have E(t — * t c , x = 0) — > 0. 
The case £ = const arises as a first integral of the motion 
in the static case and the limit £ — > corresponds to the 
singular soliton solution of Sen |16( . It is natural to ex- 
pect, then, that as t — > £ c , near x = 0, the tachyon field 
T{t = t c , x = 0) coincides with the stable kink solution of 
Sen and the time evolution in this neighbourhood stops. 
We will argue that at this point a codimension-one brane 
has formed. 



remains well-defined. Defining the momentum conju- 
gate to T as II = SS/8T the Hamiltonian is given by 
H = v/IP + V(T) VI + diTdiT. It is useful to rewrite 
Hamilton's equations of motion, n = — ^ and T — ^ , 
in a manifestly covariant way as 



d^Td^T + 1 = 



V{Tf 
S 2 



and 



E [Ed^T] = V(T)V'(T), 



(9) 



(10) 



where S = II/T 1 is defined in ©. In the limit V(T) — > 
equations ©, (0 yield 



and 



d^Td^T +1 = 



9 M [S9 M T] = 0. 



(11) 



(12) 



The solutions of (|llfl were found for arbitrary Cauchy 
data in 01 using the method of characteristics. The 
generic solutions exhibit the formation of caustics where 
second and higher order derivatives become singular. 
Caustics are known to form in systems with a pressure- 
less fluid, which is a good description of the tachyon field 
as it approaches its ground state T — > 00. It is not 
known whether the caustics are a genuine prediction of 
string theory or just an artifact of the derivative trun- 
cation which leads to the Born-Infeld Lagrangian for the 
tachyon. 

In any case, caustics are not present in the simplest 
solution of eq. i|ll|l . T = ±i, which is the asymptotic 
form for the homogeneously rolling tachyon. For T 2 = 
1, eq. I|12|l. which is equivalent to energy conservation, 
implies that £(i, x) = S(x) is an arbitrary function of 
x. In this regime the energy momentum tensor Q is 
identical to that of pressureless dust = Tj{x)u^u u 
where = d^T is interpreted as the local velocity vector 
and E is interpreted as a Lorentz-invariant matter density 



C. Vacuum Solutions 

Away from the site of the kink the field is expected 
to roll towards the vacuum T — > ±00 so that V(T) — > 
at late times for 1 / 0. To analytically study the dy- 
namics near the vacuum it is easiest to work in the 
Hamiltonian formalism [2lJ,|22j since the Lagrangian van- 
ishes in the limit V(T) — * 0, whereas the Hamiltonian 



2 For the remainder of this section the dot denotes differentiation 
with respect to time while the prime denotes differentiation with 
respect to the x coordinate. 



D. Stress-Energy Tensor 

We are interested in the behavior of Tmn as t — > t c . 
First we consider the neighbourhood near x = where 
T(t,x) = kx/(t c — t) as t — > t c (seelSJ. Near x = 0, the 
Hamiltonian is 

TV) k f k x \ 

and 

lim Too = \/naTp6(x) = t p _ 1 5(x) (13) 

t— >tc 
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using the normalization (J2J) for the potential. Similarly 



1 jlu 



Tp-1 



and Tn ^ as x — > 0. 

Consider now the late-time behavior of Tm n away from 
the site of the kink. Using the solutions of section III CI 
we find Too ~~ > 5] 0*0 while Tn, Toi and tend to zero 
for x 7^ 0. 

To summarize, we find that in the limit of condensation 
the energy momentum tensor is identical to that of a 
D(p — l)-brane: 

Too = Tp-i 5(x) + £(x) 

Tn = T i = 

T/ip = -Tp_i <5(x) <5/ii>. 

The extra bulk energy density is similar to the result 
in |l9j ] and corresponds to what has been dubbed tachyon 
matter. 



III. EFFECTIVE TACHYON FIELD THEORY 
ON THE BRANE-ANTIBRANE PAIR 

We would like to generalize the results of the previous 
section to study the dynamical formation of a tachyon 
vortex and hence a codimension-two brane. This is the 
more realistic situation, since the stable D-branes of a 
given string theory are those whose dimensions differ by 
multiples of two. We will work with an effective ac- 
tion proposed by Sen f° r the tachyon on a brane- 
antibrane pair. The field content for this system is a 

complex tachyon field T, massless gauge fields 
and scalar fields Y^, YL-. corresponding to the trans- 
verse fluctuations of the branes. The index (i) — (1), (2), 
which we call the brane index, labels which of the origi- 
nal branes (actually the brane or the antibrane) the field 
is associated with. The effective action is: 

S = - J V(T, Yfa-Yfa) (V-detMM 

+ V-detMW) d p+1 x (14) 



where 3 



M 



(0 
MN 



9mn + &'F$ N + 9mY^8 n Y^ 



1 



D M TD N T 



1 



D M T*D N T 



N-> 



(15) 



3 Greek indices {/x, u} are now understood to run over the coordi- 
nates {0, 1, 2} on which the vortex solutions depend, and hatted 
greek indices {/i, i>} run over the spatial coordinates parallel to 
the vortex {2, 3, • • ■ , p}, where p = 6 for a vortex which describes 
a 3-brane. Upper case roman indices {M, N} still run over the 
full space-time coordinates {0, 1, • ■ ■ ,p}. Finally it will be conve- 
nient later on to refer to lower case roman indices {m, n} which 
run over only the time and radial coordinates {0, 1}. 



Fmn = 9m A<# - d N A^ 



m ■ 



D 



M 



8 



M 



iA^J -(- iA^ ^ 



Ai ■ 

For the remainder of this paper we will ignore the 
transverse scalars and choose V(T, 0) = V(T) = 
t p exp (— \T\ 2 /a 2 ) where a is chosen so that the static 
singular vortex solutions of the theory Ijl4(l have the cor- 
rect tension to be interpreted as codimension 2 D-branes 
according the the normalization proposed in 11]. We 
will discuss the normalization of the potential proposed 
in [Tl| in more detail when we calculate the energy mo- 
mentum tensor for the theory (|14fl . 

Though the action l|14l) was not derived from first prin- 
ciples it obeys several necessary consistency conditions 
which are discussed in [ljj. There have been various 
other proposals for the tachyon effective action and vor- 
tex solutions on the brane- antibrane pair 23, 2J,|25j. See 
2fj] for a discussion of various models including l|14fl • 



IV. VORTEX SOLUTIONS ON THE 
BRANE-ANTIBRANE PAIR 

To construct vortex solutions we ignore the transverse 
scalars and take the remaining fields to depend only on 
the polar coordinates x^ — (x°, a; 1 , x 2 ) = (t,r,6) with 
metric g^dx^dx" = —dt 2 + dr 2 + r 2 d8 2 . Since the vortex 
solution should have azimuthal symmetry we make the 
ansatz: 

T(t,r,9) = e* s f(t,r), A$> = -A™ = lg(t,r) (17) 

with all other components of A^ vanishing. This gener- 
alizes the ansatz used in 1 1 if to include time-dependence 
in the fields / and g. For (|17(l one has 

D t T = e w f(t,r), D r T = e ie f'(t,r), 



D e T = e ie z(l-g(t,r))f(t,r) 



and 



F, 



- „ F ( 2 ) 



where the dot denotes differentiation with respect to time 
and the prime now denotes differentiation with respect 

to the radial coordinate. The matrices M^ N are 



M 



(i) 

MN 



M 



(2) 
MN 



-1 + / 2 //' a'g/2 

//' 1 + f 2 a'g'/2 

-a'g/2 -a'g'/2 r 2 + (1 - g) 2 f 2 

6^ 

T 



M 



(i) 

MN 



(18) 



We also have det(AfW) = det(M< 2 )) since M$ N = 

MjyL and so we omit the brane index (i) on det(AfW) 
in subsequent calculations. 
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The action for this ansatz simplifies to: 



V(f)((l + d m fd m f) [r 2 + / 2 (l- 5 ) 2 



/2 12 \ 

+ ^d m gd m g-^(e mn d m fd n g) 2 ) <F+^19) 

where x rn = (x^^x 1 ) = (t,r) and g mn dx m dx n = —dt 2 + 
dr 2 . For notational convenience we define the scalar 
quantity 

V(f) 



S(t,r) 



(20) 



yf- det(M) 

in analogy with jnj). The equation of motion for the 
tachyon is 



= £(l + d"7cW)(l-.g) 2 / + 



^2 



(21) 



and the nontrivial component of the equation of motion 
for the gauge field is 

^-d m [m m g ~ T,{e ab djd b g)e mn d n f] 

= E(l + d m fd m f)f(g-l). (22) 

Alhough these equations are somewhat cumbersome, in- 
spection of 1221) tells us that there should exist a solu- 
tion g(t,r) such that at g = 1, the vacuum, we have 
d m g — 0. This is the asymptotic behavior which cor- 
responds to a vortex solution; it is already known from 
that the static solution g(r) is a monotonically in- 
creasing function which varies between and 1. Thus we 
shall only consider solutions with the asymptotic behav- 
ior d m g(t,r) -40 as g{t 7 r) — ► 1. We will take initial 
data g(0, r) = <?j(r) such that <7i(0) = 0, < gi{r) < 1 for 
all r and g[{r) > for all r. In addition we will focus on 
initial tachyon profiles /(0, r) = fi(r) such that fi(Q) = 
and fl(r) > for all r. For these initial conditions the 
tachyon must start rolling for r ^ due to its displace- 
ment from the unstable vacuum V'(f) = 0. Since the 
asymptotic gi{r — > oo) — ► 1 is an exact solution of 1221) 
and g(t, r = 0) = by construction, we therefore expect 
g(r, t) to increase towards unity for finite r ^ 0. 

A. Solutions Near the Core of the Defect 

As in ^3 > to analytically study the dynamics near the 
core of the vortex, r = 0, we make the ansatz: 



f{t,r)^p(t)r, g{t,r)^q(t)r 



(23) 



for small r. Dropping terms which are sublcading in r 
yields a set of coupled ODEs for p(t) and q(t) 



pq 4 q—pq' l +2pq i q—2pq 2 q 2 + 



T p 3 q 2 - 



tPQ 



tPQ" = 
(24) 



6|— ' 1 ' 

5- 

CT< - I p(t) 

4 - ■-■ q(t) 
0, - 1/P(0 



\ ' 
\ 



3- 



l/q(t) 



2- 



0.01 0.02 0.03 0.04 0.05 0.06 
t 

FIG. 1: Numerical solution for p(t) and q(t) of 1241251 1 show- 
ing the finite-time slope divergence. 1/p and 1/q are also 
shown, demonstrating the linearity of these functions near 
the critical time. 



from and 



•• 3 i 2 2 •• , 2- n -2 , o • 2 • n 2 -2 , ^ 4 
ppq + pqq + qq-2qq + 2ppq q - 2p qq + —p q 



+ -ptp q + -~q + -sp r = o 



(25) 



from (J22J). Although I|24I25|I are difficult to solve ana- 
lytically it is straightforward to verify that in the regime 
where p, q and higher derivatives are large compared to p, 
q there exists an approximate solution to I|24I25|I where 
both p(t) and q{t) are divergent in finite time t c : 



Pit) 



Po 

tr-t' 



?(*) = 



</<> 



tr-t 



(26) 



in analogy with the kink solution of |Ioj|. Numerical so- 
lutions to (|24I25() agree with this prediction, as shown in 
figure 



B. Solutions Away From the Core of the Defect 

We are interested in solutions where d m g — ► as g — > 1 
and where f(t, r) — * oo as t — ► oo. Since, as we have seen 
above, g'(t,r — 0) is diverging in finite time, therefore 
g must be increasing to unity for r =/= so that at late 
times we expect g(t, r) to resemble a step function. Thus 
to study the dynamics away from the core of the defect 
we begin with the ansatz: 



g (t,r) = 1 - ea(t,r) 



(27) 



and work only to leading order in e. The leading order 
contribution to (|21|) decouples completely from er(t, r) 



7(i + n /"(i - h vf? %(i -i 2 + n 



i-/ 2 + /' 2 J' = o 



(28) 
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We can consistently find solutions of l(2*5)l by taking / to 
be a solution of the eikonal equation 1 — f 2 + f' 2 = 0. 
Subject to this constraint the second term in the square 
braces in (|28|l vanishes trivially. The constraint that the 
first term in the square braces in 1|28J) vanishes is exactly 
the same as the equation of motion one would derive from 
C = —V{f)\/1 + d m fd m f, the Born-Infeld Lagrangian. 
The eikonal equation yields the Born-Infcld equation as 
a differential consequence, which is not surprising since 
this amounts to minimizing the action by setting C = 0. 
Thus the PDE l|28l) is automatically satisfied when / is 
a solution of the eikonal equation. We find, then, that 
when g(t, r) = 1 the tachyon field must obey 



d m fd m f + 1 = 



(29) 



as in section ITT1 

The ansatz (|27|) yields no simplification of 1)22(1: how- 
ever, we may solve for a(t, r) given l|29|l based on more 
fundamental constraints. The argument of the square 
root in l|19|l must be nonnegative to ensure reality of the 
Lagrangian. Thus for f(t, r) given by i|29|) the require- 
ment that the Lagrangian be real translates to 

d m ad m (7 - (e mn d m fd n a) 2 > 0. 

or, with / given by (29). 



cT^l + f l2 -(j'f) >0. 



(30) 



For real fields it is clear that 1|30[) can only be solved when 
the equality is taken. It is worth noting that since the 
Lagrangian vanishes when the equality is taken in (|30|l . 
this constraint ensures that the full equations of motion 
are satisfied. 

As in section |n] we can avoid the difficulties of caustic 
formation in the general solutions of (|29|l found in 01 
by taking the Cauchy data to be linear and using the 
one-parameter family of solutions 



f(t,r) =at+ Va 2 - I r. 
Reality of the Lagrangian requires 



(31) 



la' = 0. 



This PDE is separable and we find the solution 



g(t,r) 



e exp 



R 



1 



(32) 



where R is a separation constant. Note that the solution 
(|32|l becomes static in the limit a 2 — > 1, the homogeneous 
rolling tachyon. In fact, for a 2 = 1 any function a(t, r) = 
cr(r) satisfying the necessary boundary condition a(r — > 
oo) — > will generate a solution. We will ultimately be 
interested in this limit. 4 



It is noteworthy that taking g(t, r) close to unity 1271) 
ultimately translates into the requirement that det(M) 
must vanish. The solutions I(31I32|I should be thought of 
as late-time asymptotics where V(f) — > since / — > oo. 
In this limit £(f , r) defined in 120|) has the indeterminant 
form 2 as in section^ The quantity E is of some interest 
for two reasons. First, it parametrizes the manner in 
which we take the limits V(f) — > and det(^4) — > as 
we approach the vacuum state. Second, the form of E for 
r =/= will determine the form of the energy-momentum 
tensor in that regime, as in the case of the kink. Following 
the discussion in section[n]we will use energy-momentum 
conservation to place constraints on the aymptotic form 
of E. 



V. STRESS-ENERGY TENSOR 

In this section we demonstrate that the vortex solu- 
tions found above give rise to the formation of a sin- 
gularity in the stress-energy tensor of the tachyon field, 
which corresponds exactly with that of a codimension- 
two D-brane in the final state, whose tension has the 
value expected from string theory. We also derive the 
bulk stress-energy tensor for the leftover tachyon matter, 
which continues to roll even after the formation of the 
D-brane. 

The stress-energy tensor for the action (|f4H is 



T 



MN 



■(2)1 



V / -detM( 1 )(M ( " ) 1 ) 



MN 



+ v / -dctM( 2 )(M ( - ) 1 ) 



MN 
(2)) S 



(33) 



where the subscript S denotes the symmetric part of 
the matrix, i.e., (M^)f N = \[{M~^) MN + (M^) MN }. 

The components of T MN parallel to the vortex simplify 
to 



'-V(f)y/~det(M)5^. 



(34) 



For the components involving i, r, 0, it is useful to rewrite 
T^ u in terms of E and the symmetrized cofactor matrix 
of M$, which we define as C^ 5 Since Cfii = C%. we 
drop the brane index on the cofactor matrices. In the 
t, r, 6 directions, then, we have 



2E 
r 



(35) 



The nonzero components of the cofactor matrix for the 
Lagrangian l|19() are 



C tt = [r 2 + f{l- g ) 2 }{l + f' 2 ) + l - T g' 2 , 



4 The exact functional form of a at large r and late times turns 
out to be of little importance to the ensuing analysis. 



5 That is to say C, M T 



:det(Af)(M- 1 )^. 
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-[r 2 + f(i-g) 2 ] ff'-^-gg', 
-V 2 + f 2 {i- 9 ) 2 ] (i-/ 2 ) + ^<? 2 , 



c ee = _ (1 _ /2 + / / 2) 



A. Normalization of the Potential 

In ^3 Sen finds that the action l(Ti|) provides a good ef- 
fective description of the tachyon on the brane-antibrane 
system provided the potential V(T) is chosen to satisfy 
the normalization constraint 



Tp-2 



= -It / V(z)^z 2 {l - G{z)) 2 + — (G'{z)) 2 dz 

(36) 

where r p _2 = (27r) 2 a'r p is the tension of a (p — 2)-brane, 
G(z) = G(F~ l (z)), and {/(r) = F(br), g{r) = G(br)} 
are the static soliton solutions to be understood in the 
limit that b — > oo. This constraint is necessary to ensure 
that the vortex solution has the correct tension to be 
interpreted as aD(p- 2)-brane. In the time-dependent 
case there is some ambiguity as to how to interpret <|36[) 
since this statement appears to depend on the functional 
form of the solutions, which would make the right-hand- 
side apparently time-dependent. 6 

But physically, it makes sense to impose <|36|) at t > t c 
since t c is the time by which the brane has actually 
formed, and in the limit t — > t c the time-dependent so- 
lutions should coincide with the soliton solutions in the 
neighbourhood of r = 0. In the time-dependent case, it is 
(t c — t)^ 1 which tends to infinity (as t — > t c ) and plays the 
role of b. Although the exact functional forms of G(z), 
F(z) are not known for all z, we can infer from (|26[l that 
G{z) = qoz and F(z) = poz near z = 0. Furthermore, we 
know that G(z) — > 1 for sufficiently large z by construc- 
tion. We also have F(z = 0) = and F(z ^ 0) ^ so 
that F-^O) = and F'^z ^0)^0. 

Let us consider the two terms under the square root 
in 136|) . The first term, z 2 (l — G(z)) 2 , is small near 
z = due to the overall multiplicative factor of z 2 . On 
the other hand, at large z this term is also small since 
F~ x {z) is large and hence G{z) = G(F~ x {z)) = 1. We 
conclude that the derivative term under the root in i|36|) 
dominates. At small z we have F^ 1 (z) = p^ 1 z and thus 
G(z) S qoF^iz) S 3±z so that G'(z) = ^. For simplic- 
ity we take G'(z) = ^ for all z since this expression is 
multiplied by V(z) which tends to zero quickly for large 
z. We find then that 



So 



Pa J- 



2 / 2 

" z l a dz = 4tt. 



(37) 



Arguments are presented in ll ll for why 1361 is in fact only a 
constraint on V{T) and not on the solutions T(x) themselves. 



The normalization l|37() is equivalent to a = 4y / 7r|^. We 
shall see later on that the relation a = 4y / 7r^- may be 
equivalently viewed as a constraint on the arbitrary func- 
tion E(t, r). 



B. Stress-Energy Tensor at r — 

At r = and t — > t c the solutions Ij26(l are valid and 
the Hamiltonian is 



rpU 



,/2 



^([r 2 + / 2 (i-. 9 ) 2 ] d + n + ^g' 1 



Tp gg a 
r t r — t 



exp 



2 2 



to leading order in r. Then, using a = 4y / 7r^j-, 



lim T u 



S(r) 



— Tp-2 S(r cos 6) 8{r sin 9) 
and the components parallel to the vortex are 



--S^V(f)^-det(M) 
-t p -2 S^Sir cos 6)6(r sin B). 



The remaining components, T 11 and T 01 , are vanishing 
at r = 0. The angular component T 22 = T ee contains 
a delta function at r — 0; however this is an artifact of 
6 being a bad coordinate at r = 0; it can be seen that 
Tge = 0, and going to Cartesian coordinates confirms 
that the T^ u = for the transverse coordinates, as should 
be the case for a D-brane. This result is the same as in 
the static case 10. 



C. Stress-Energy Tensor at r > 



For r > at late times the solutions (131132(1 are valid. 
For simplicity we take a 2 = 1 and work only to leading 
order in e. The Hamiltonian is 



T 



2r£(t,r) 



and the remaining components of T MN vanish at late 
times for r > 0. Notice that conservation of energy 
q mT mn = o at large r forces 

E(t,r)=E(r). 
That is, £ is an arbitrary function of r as in section ITT1 

To summarize, we find that in the limit of condensation 
the energy momentum tensor is identical to that of a 
T)(p — 2) -brane 



T 00 = T p - 2 6(r cos 6) 6(rsm8) + 2r£(r) 

/yill rj-i'2'2 Q 

T"* = -5> iP Tp-26(rcos6)6(rsm6) 



(38) 
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with all off-diagonal components vanishing. The extra 
bulk energy density 2r£(r) is similar to the result in [l^ 
and section |n] and corresponds to tachyon matter rolling 
toward T — ► oo in the bulk. 



D. Conservation of Energy 

We can constrain E(r) using conservation of energy. 
Initially the system consists of two Dp-branes with energy 
density 2t p Vi in the (p— 2)-dimensional subspace spanned 
by {x^ 1 }, where V2 is the volume of the 2-dimensional sub- 
space spanned by {r, 6}. At late times, after the codimen- 
sion 2 brane and its antibrane have formed, the energy 
density in the (p — 2)-dimensional space is given by the 
sum of the D(p— 2)-brane tensions, 2r p _2, and the energy 
density due to tachyon matter A-k J r 2 E(r) dr. Conserva- 
tion of energy thus implies 

2t p V 2 = 2t p _2+47t J drr 2 Y,(r) 

(2F 2 - 2(2tt)V) t p = 4tt J drr 2 E(r). 



The matrix M$ N of eq. JTJSJ can be written in block 



Since £(r) is arbitrary we can take 



r 4irr 2 



where S(r) satisfies the constraint 
dr S(r) = 



■K 2 CtTr> 



(39) 



(40) 



The conditions l(39|l . 1(40(1 are equivalent to 1(37(1 and may 
be thought of as an alternative to the normalization ((36(1 . 



VI. INCLUSION OF MASSLESS GAUGE 
FIELDS 

We will now restrict ourselves to a (5+l)-dimensional 
spacetime with {M, N}={0, 1, ■ ■ ■ , 5}, {fi, v}={0, 1, 2} 
and {/}, £}={3,4, 5}. There are two gauge fields in the 
problem: A¥L and A%L , or equivalently A^ 1 = A%L 

and At 1 = A^ — A^L, which have different couplings 
to the tachyon. We have already shown that At is the 
field which condenses in the vortex, hence its associated 
gauge symmetry is spontaneously broken. For reheating 
it is thus A+ which most closely resembles the Standard 
Model photon. We will ignore fluctuations of the heavy 
fields At, A+, and At, keeping only the background so- 
lution for At (which was given in section IIV|) , and the 

fluctuations of the photon At. This leads to considerable 
simplification since it ensures that D^T — 0. To compute 
the production of photons in the time-dependent back- 
ground, we want to expand the action (|14|l to quadratic 
order in At- 



diagonal form as 



M 



(i) 
MN 



V, 



« 



s 



« 



s 



fJ,V 



(41) 



where V^> = M$, is the contribution from the vortex 



background given in (|18|l . S^l — a'd^A^ 1 is the contri 



bution from {t, r, 9} derivatives of A^\ and F¥J is the 

field strength tensor for A^ . Using a well-known identity 
for determinants we can write 

det(M^) = det V (l) det (l + afF® + S^V^S^ 

(42) 

Expanding det(l + a'F^ + S^V^S^) to quadratic or- 
der in A^ , the action (|14(l becomes 



S S 



dpA+d^ 



d p+1 xV(f)V^drtV (iV-^d^AtduA* 

(43) 



where we have chosen the gauge dp, At = and disre- 
garded the piece which does not depend on A+ . We omit 
the brane index on Vn) since the determinant and the 
symmetric part of V^ 1 are equal for both i = {1,2}. 
Defining an effective metric Gmn by 



[G 



MN\ 



V (S) 

' [0] 



[0] 



(44) 



the action ((43(1 may be written as 




ac MN S^d M Atd N Af dP +1 x. 

(45) 

From 1(45(1 one sees that the fluctuations of the photon 
behave like a collection of massless scalar fields propa- 
gating in a nonflat spacetime described by the metric 
Gmn, with a position- and time-dependent gauge cou- 
pling given by g 2 = l/V(f(t, r)). 

To get an intuitive sense for the behavior of the ac- 
tion (|4"5|) we note that the stress-energy tensor derived in 
section can be written as 



T 



MN 



-v(fh 



-iMN 



In the limit of condensation, T MN is given by l(38|l . so 
that once the brane has formed the action 145(1 reduces 
to a description of gauge fields propagating in a (3+1)- 
dimensional Minkowski space, with an additional compo- 
nent which couples the gauge fields to the tachyon matter 
density in the bulk. In other words, the effective metric 
Gmn starts off being smooth throughout the bulk, but 
within the time i c , its support collapses to become a delta 
function S^ 2 \x) in the relevant extra dimensions {r, 8}. 
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The equations of motion resulting from the effective 
action (|45|l are difficult to solve analytically since the 
effective metric Gmn depends nontrivially on both r and 
t and is nondiagonal in the subspace of {t, r}. For this 
reason we would like to propose a simplified model of 
the condensation which captures the essential features of 
the action 1|45[) . We have derived solutions for the vortex 



background valid at small r, 



(t c — i), and at large r, 



t ^ (t c — t). Similarly, the energy momentum tensor we 
have derived corresponding to these solutions has very 
different behavior in the r < (t c — t) and r > (t c — t) 
regions of the spacetime. 

For r < (t c — t) the energy momentum tensor contracts 
to a delta function centered at r — 0, with (t c — t) play- 
ing the role of the small parameter which regularizes the 
delta function. That is to say, 



tT mn ^ 



1 



tr.-t 



exp 



2 2 



a 2 (t c - ty 



H 



MN 



Finally we argue that the bulk part of the action can 
be ignored. To this end, let us change to coordinates 
which are comoving with the contraction of the vortex 
core: 

r 



tr-t 



t = t r 



t. 



(46) 



In terms of these coordinates the "small r" solutions are 
valid for f < 1 and the "large r" solutions are valid for 
f > I. The Jacobian of this transformation is —t so that 



Sbranc OC 

exp 



2 ~2 



H MN S^d M Atd N At dtdrd6dx & . 



To lowest order in f 2 the matrix entries H MN in terms of 
these new coordinates are all constant. Consider now the 
piece of the action which couples to the tachyon matter 
density in the bulk, written in terms of these new coor- 
dinates: 



at small r, where the matrix entries H 00 , are finite 
as t — > t c and the remaining components of H MN tend 
to zero (near r = 0) as t — > t c . 

In the r > (t c — t) region, the energy momentum tensor 
has quite different behavior. After condensation of the 
defect has completed at t = t c , the energy density in the 
bulk (r > 0) is due entirely to tachyon matter, while the 
part of the stress-energy which is going into the tension 
of the defect vanishes in this region. Hence 



rT 



MN 



2«r 2 £(r) 



as t — > t c , for r > (t c — t). 

In our simplified model of the particle production 
due to the tachyon condensation we therefore split the 
energy momentum tensor into brane and bulk pieces 
T MN = T«~ c + where T-~ c contracts to a delta 

function as t -> t c and T b A ^ -> 26tf6$rE(r) in the same 
limit. The action (|45|l then splits into two components 
S = 5buik + <Sbranc- We expect most of the particle pro- 
duction to occur near the end of the condensation, when 
the background tachyon field is becoming singular near 
the vortex, so the best approximations for the simpli- 
fied gauge field action are those which describe the exact 
expression most accurately near t = t c : 



Sbrane ^ 



1 



t r - t 



exp 



2 2 



a 2 {t c - ty 



H MN 5 p ' i> dMA±d N A± dtdrd9dx & 



and 



Sbulk oc 



(r) (- 



dtdrd8dx a . 



Si 



bulk OC 



PrE(ft) -5^ At At dtdrdOdx 



Since E(z) is not singular at z = 0, the bulk piece of the 
action become negligible near the end of the contraction 
t — > 0. This is a consequence of the fact that as the 
condensation proceeds the gauge field is confined to the 
descendant brane. 



VII. SIMPLIFIED MODEL OF REHEATING 

In the previous section we argued that the gauge field 
couples most strongly to the part of the tachyon back- 
ground which is collapsing to form the defect. This 
closely resembles a gauge theory defined on a manifold 
in which a two-dimensional subspace which is shrinking 
with time. As a simplified model of the interaction we 
thus consider a massless spin- 1 field 



~ MA g NB F MN F AB d? +1 x 



propagating in a FRW-like background 



gMNdx M dx N 



-dt 2 + R(t) 2 (df 2 +f 2 d9 2 )+c 



b( L ydx ll dx v . 

The coordinate f in (|47|) is fixed with the expansion and 
is thus corresponds to f defined in (|46|) . However for 
simplicity of notation we will drop the tilde and write r 
instead of f in the remainder of the paper. We take r to 
be dimensionless while t, x^ and R have dimensions of 
length. 

If we restrict ourselves to configurations with A^ = 0, 
Ap- ^ o 7 and impose the gauge condition d^A^ — then 



At earlier times the coefficient of the bulk part of the La- 
grangian would have time dependence, and the bulk La- 
grangian would contain contributions from all the deriva- 
tive of the gauge field, but this form is valid close to t c . 



7 This restriction will only underestimate the reheating. Since we 
want to show that the reheating can be efficient, this approxi- 
mation will not weaken the ensuing argument. 
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the action simplifies to 



The particular solutions of ifl^l are labeled by the mo- 



-4 



'99 



„MN 



menta k in the 3 large dimensions ({x , a 
^dMAjxdxAp d p+1 x (48) take them to be normalized according to 



'}) and we 



Notice that for the metric (|47(l and the gauge field config- 



uration = we have chosen, one has = 9m Au 



and R MN A M A N 



RmAi>A u = 0. 



We will impose homogeneous boundary conditions at 
r = 1 and take the scale factor in 147fl to be 



R(t) 




if t <0; 
if < t < t c 
if t > t r . 



(49) 



where Rq represents the initial radial size of the extra di- 
mensions. This approximation for the time-dependence 
of the vortex core is the simplest form which has the same 
qualitative behavior as the true background, while still 
allowing us to solve analytically for the gauge field wave 
functions in the background. A shortcoming of this ap- 
proximation is that R is discontinuous at the interfaces, 
which leads to an ultraviolet divergence in the produc- 
tion of gauge bosons. The behavior of the actual R(t) 
is smooth, and must yield a finite amount of particle 
production |{|. In addition, we will find a separate UV 
divergence in the particle production in the limit as the 
vortex core thickness goes to zero. This is presumably an 
artifact of the effective field theory which is not present 
in the full string theory, and we deal with it by intro- 
ducing the cutoff e on the final radius of the defect core, 
which we will take to be of order the string length l a . 



A. Gauge Field Solutions 

The first step in computing the production of photons 
in the time-dependent background is to solve their equa- 
tion of motion following from 



.1? 



Ap" Ap' dlAP 



- 2 R AP - AP+ iv + ^ + i^ +d ^ AP = (50) 

where the dot and prime denote differentiation with re- 
spect to t and r, respectively. 
Equation (|50|l separates as 



AP{t,r,6,x*) = <j>{t)<fi(r)e{6) x P(x fl ) (51) 



where 



d%x" = -*V, dox° = 0, k 2 = k^k^ = k-k (52) 

c^6 = -m 2 e (53) 

^ + V+(c 2 -^ = (54) 

^ + 2 ^+ (^ + fc2 )^ = ° (55) 
and c is a separation constant. 



hp- 



The particular solution of (|53|l is a sum of sines and 
cosines. The odd parity and even parity modes (under 
6 — > —8) do not mix with the even ones, and for simplic- 
ity, we restrict our attention to the even modes, which 
include the massless one. This can underestimate the ef- 
ficiency of the reheating by a factor of 2 at most. The 
solution of (|53l) is thus 



e m (0) 



i 



cos(mO). 



These are orthogonal for different values of m, and re- 
quiring that the solution be single-valued restricts m to 
be an integer. 

So far we have not been specific about the geometry 
of the two extra dimensions around which the original 
annihilating branes were wrapped. One simple possibil- 
ity is a 2-sphere, where the descendant brane and an- 
tibrane (vortices) form at antipodal points. Since all the 
singular behavior of the tachyon background is localized 
near these points, the curvature and topology in the bulk 
should have little effect on particle production near the 
defects. To simplify the mathematics, we therefore re- 
place either of the two hemispheres of the sphere with 
a flat disk, in the coordinate region r < 1. The correct 
boundary condition on radial eigenfunctions of the bulk 
Laplacian is that their derivatives vanish at r = 1, so that 
they are smooth at the interface where the two halves of 
the space are glued together. 

The solution of subject to the boundary condition 
tp'(r = 1) = and the requirement that (p (r) be regular 
at the origin, is 



V2 



Jm-\-l (Cmn) 



J rn. (CmiJ 1 ) 



(56) 



where c mn is the nth zero of J' m (r). The solutions l|56l) are 
orthogonal for different values of n. The zero mode n = 1, 
m = must be treated separately; it is the constant 
solution, where Coi = and </ 3 oi( r ) = 1- 

The solution of (|55|l depends on the scale factor. For 
t < and t > t c the solutions are trivial and are given 
by 



(/>mn{t) 



1 



RoV^ 



and 



e\/2ui„ 



respectively. We have defined u 



2 

inn 



•*) (57) 



(58) 



+ k 2 and 



k 2 . The multiplicative factors in (|57|l and 
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(|58|l are introduced for later convenience and ensure that 
a mn and d mn will be properly normalized annihilation 
operators in the appropriate spacetime region when the 
gauge field is quantized. It will be convenient for what 
follows to introduce phase-shifted annihilation and cre- 
ation operators in the region t > t c : d m n — e~ l ^ mnta d m n 



and c4 n = 
becomes 

(t>mn{t) = 



1 



€y/2uj mn 



c dJ 7ln . In terms of these operators l(58j) 



l%m c ' "(Tin I 



We are suppressing the dependence of the annihila- 
tion/creation operators on the 3-momenta k. 

In the region < t < t c where the scale factor depends 
nontrivially on time, the solution of (|55|l is 



4^mn (^) 



H - Crnn J —z 



'■{Ro-Tft) 



'■[Ra-rft) 



(59) 



where p mn = ^- \/r) 2 — 4c^„ . Some comments are in or- 
der concerning this solution, which has different behavior 
for massive and massless modes. 



In the massless case c mn — 0, p mn — 1/2 and J p , 
J_ p are linearly independent, since p is noninteger. The 
constants B mn and C mn in fflfy should be interpreted as 
independent, real- valued constants so that 4> mn (t) is real. 

In the massive case c mn ^ 0, rj 2 < l, 8 and the or- 
der of the Bessel functions in l(5*9*jl. p mn , is pure imagi- 
nary. Such Bessel functions are complex- valued and there 
are several options for constructing real solutions. Since 
(J v (x))* = J v * (x) for real x, we can impose C mn = B} nrl . 
Calculating the Wronskian of Ji S , J-i S verifies that for 
real s ^ these two solutions are linearly independent. 
Since we do not need to explicitly quantize the field in 
the region < t < t c , the dagger can be thought of sim- 
ply as complex conjugation. The interpretation of B mn 
and B\ nn as annihilation and creation operators is un- 
necessary, since there are no asymptotic states in this 
region. 

We summarize this subsection by putting these results 
together to write the general solution of 150(1 as 

AP(X M ) = £Hx§(z A )em(0)0m„W¥>mn(r). 



where 



1 I A p-«S mn (t-t c ) _|_ Jt p- 



(t-tc)) 



if t < 0; 
„[^(Ro-vt)}) if0<<<i c 
if t > t c . 



B. Spectra of Produced Particles 



Massive Modes 



The next step is to impose continuity of and 8m A^ 
at t = and t = t c in order to compute the Bogoliubov 
coefficients, which relate the annihilation and creation 
operators for t > t c (d mn and efj^) to those in the region 
t < (a mn and a}„„). 9 Smoothness of the solutions 
at the interfaces is ensured by the continuity of <j) mn and 
4>mn, since 4> mn is the only part of A p which changes 
between the different spacetime regions. 



One might expect rj 1 < 1 on physical grounds. In fact, for the 
ensuing arguments to hold one need only restrict r\ < 3.7 to 
ensure that p m n is pure imaginary and nonzero for all massive 
modes. 

We are ultimately interested in calculating the number of d mn 
quanta in the vacuum annihilated by a mn . For this purpose it is 
just as good to use d mn , dln n as d mn , d} mn since the phase shift 

mn u, mrt 



cancels out of the number operator, d. 



For the massive modes continuity of </> mn (i) and </> mn (i) 
at t = implies 



= U 



Ml U 2 

,t 



Un U\ 



it ~ \ 

mn / \ ran 

The entries of the matrix U are given by 

ir] 



fit 

mn 



(60) 



Ul = 



' R U! ri 



U 2 = 



\\ 2 2y/2R~^r ri 
ik x f^J' Pmn (kRo/rj) 

V ^^mn 

irj 



2 v / 2i?o ti; ri 



- ihf^-J'- Pmn (kRo/ V ). 

V ^mn 

Continuity at t = t c similarly gives 



III II 

"t 

III II 



V 



B m 
Bt 



Ul v 2 



J Pmn {kR /r)) 



J- Pm JkRo/r]) 



-E$m 



(61) 



12 



The matrix entries Vi are obtained from m by replacing 
i?o with e and uj mn with u) mn - 
From 1)60(1 . I(61|) we can write 



= Vil- 



li 



The last equality defines the Bogoliubov coefficients. 
Note that there is no summation implied over any of 
the indices in the above expression. The indices m and n 
label the modes of the in- and out-states, and the sum- 
mation which appears in the general definition of the Bo- 
goliubov coefficients is not present here. 

We can now determine the spectrum of Kaluza-Klein 
(KK) excitations of the photon which is produced in the 
tachyon vortex background. Observers in the future see 
a spectrum of massive particles in the final state given 
by 



I A, 



1 



\detU\ 



:\uiV2 - U2V1Y 



The determinant is 



det U = detV = -sin(7rp mn ) 



(62) 



which may be obtained by using the Wronskian of J p and 
J^ p . The fact that the two determinants are equal en- 
sures the appropriate normalization \a mn \ 2 — \/3 m n\ 2 = 1 
of the Bogoliubov coefficients. The explicit expression 
for iV^™ (fc) can be obtained analytically but it is com- 
plicated and we do not write it out here; instead we will 
give numerical results. 

The mass of a KK mode with quantum numbers m, n 
is c mn /e, which increases with m and n. Figures [20 il- 
lustrate the dependence of Nffi£ Q (k) on the 3- momentum 
k for the lightest few massive modes. The parameters of 
the vortex background are taken to be Rq — 10 l s , e = l s 
and 77 = 1. (The dependence of N^ I 7 I >0 (k) on Ro, e and 77 
is essentially the same as that of the spectrum for mass- 
less modes, NM=a(k), which we will discuss in the next 
section.) 



2. Massless Modes 

The analysis for the massless modes proceeds similarly 
to the calculations in the preceding section, but more 
simply since in this case p mn = 1/2 and uj mn = u mn — k. 
The Bogoliubov coefficients do not depend on the mode 
indices: 
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The observed spectrum of particles in the final state 
N M =o(k) = \(3\ 2 (63) 
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FIG. 2: A r J2>o(' c ) versus k for increasing values of the radial 
quantum number n showing decreased production of heavier 
modes. 
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FIG. 3: Same as fig. [5] but varying the angular quantum 
number m. 



where the second line shows the leading small-e behavior. 
The exact result can also be obtained in closed form. 
Using explicit representations of J n /%, we will be able 
to integrate this expression exactly to obtain the energy 
density of produced radiation. Figures 1 4161 show plots 
of NM=o(k) as a function of k for various values of the 
parameters Ro, e and 77. 

C. Energy Density 

To find the total energy density of radiation produced 
by the massive modes we should integrate over all k and 
sum over the mode indices m and n and the number of 
polarizations (3 for massive vector bosons): 



Pm>o — 3 



d 3 k 
(2tt) 3 



k 2 + ^N^l (k). 



(64) 
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FIG. 5: A r M=o(fe) versus k for different values of Rq. e — l 3 
and T) = 1. 
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For the massless modes we only integrate over k and sum 
over the 2 polarizations, since NM>o{k) has no depen- 
dence on the mode indices: 

/d 3 k 
■p^kN M=0 (k). (65) 

As anticipated above, the integrals in (16411 . (I65|) are not 

convergent due to the fact that the simplified tachyon 
background (|49l) has discontinuous time derivatives at 
t = and t = t c . This has the consequence that the 
spectrum N(k) decreases like k~ 2 for large k, which is 
too slow for convergence. In reality nth derivatives of 
R(t) do not exceed 0(1//™), so the production of modes 
with k > I J 1 should be exponentially suppressed. We 
therefore introduce a UV cutoff, k max — A ~ Zj 1 . More- 
over N(k) is divergent as the vortex radius e — * 0, so e is 
also presumably limited by the string scale, e ~ l s . 

In the final state, the heavy KK modes have mass given 
by c mn /e, where the c m „'s are order unity and larger. 
These are near the cutoff, so their contributions are of 
the same order as other UV contributions which we are 
omitting. For consistency we should thus neglect the 
massive states' contribution to the total energy density 
on the vortex. Again, this underestimates the efficiency 
of particle production and makes our estimates conser- 
vative. Henceforth we will refer to pm=q as simply p. 

Although the integral (|65|l can be performed ana- 
lytically, the resulting expression for p is cumbersome. 
Rather than write it out explicitly we will discuss some 
noteworthy features, plot p with respect to the parame- 
ters of the model and present some useful simplifications 
of the complete expression in various limits. 

Figures !7l5l show the dependence of p on the initial size 
of the brane. For Rq greater than a few times l s , the en- 
ergy density is relatively insensitive to changes in Rq. We 
can therefore reduce the dimensionality of the parameter 
space by simply assuming that the extra dimensions are 
somewhat larger than l s . In fact in the limit of large Rq, 
the energy density takes the very simple form 

r] 2 A 2 

llm P ~ 0Z2 ( 66 ) 



which is the main result of this section. We recall that 77, 
which parametrizes the speed at which the vortex forms, 
is predicted from eq. (|46|l to be 77 = 1. 

As a check on our calculations, we have also consid- 
ered the limit as Rq — > e, which corresponds to a static 
background, with no vortex condensation. As expected, 
the energy density of produced particles goes to zero, 

A 4 , .0 A 4 
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D. Efficiency of Reheating 



FIG. 6: A r M=o(fe) versus k for different values of r\. Rq — 10 l s 
and e = l s . 



To quantify the efficiency of the reheating we need to 
determine how much energy is available to produce the 
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FIG. 7: p versus compactification radius Ro for different 
values of e = A -1 with r) = 1. Dotted line is the critical 
density for efficient reheating p c , l|67jl . 
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FIG. 8: Same as figure |7| but for different values of vortex 
collapse rate r\ and with e = A~ = 0.7 l a . 



photons on the final-state 3-brane. Initially the system 
consisted of a D5-brane plus antibrane, whose 3D en- 
ergy density was given by 2T5V2, where T5 is the ten- 
sion of a D5-brane and V% is the volume of the compact 
2-space {r, 8} wrapped by the branes. The final state 
consists of a D3-brane/antibrane with total tension 2r 3 . 
Until now we considered just half of the 2-sphere and fo- 
cused on a single vortex located at r = 0. Conservation 
of Ramond-Ramond charge requires the second vortex, 
which we place at the south pole of the sphere to preserve 
azimuthal symmetry. These two defects are identical and 
are matched at the equator of the sphere. The vortex at 
the south pole represents the D3-antibrane. 

Since reheating on each final-state brane should be 
equally efficient, the 3D energy density available for re- 
heating on one them, which we call the critical energy 
density p c , is half the difference between the initial and 



final tensions of the branes and antibranes: 
Pc = T5V2 - r 3 = r 3 I . . - 1 



V 4tt 2 



a' 



where we have used the recursion relation 2tt\ 



_i. The tension of a D3- brane is given by |27 
1 1 



73 



g s (27r) 3 a /2 ' 



The string coupling g s , in 5+1 dimensions of which two 
are compact, is determined by the gauge coupling evalu- 
ated at the string scale, a(M s ) 0: 



Vn 



27T 2 0! 



7 a(M s ). 



Thus we have 

l 4 sPc = 



167r 3 a(M s ) 



1 - 



4^ 2 



V 2 



(67) 



In the regime where V2 ~ 2ttRq ^> l 2 s , which was where 
we could most easily quantify the particle production, 
the second term in parentheses can be neglected, and 
in any case it would be unimportant for a rough esti- 
mate unless it accidentally canceled the first term (1) 
to high accuracy. Hence we drop this term and take 
l*p e = (W7r s a(M s ))-\ 

The energy density p c is the critical value at which the 
conversion into radiation would be 100% efficient. In our 
analysis we take a(M s ) = [2( which gives p c l A s = 0.05. 
The critical energy density is shown as a dashed horizon- 
tal line in figures [7| and [H] We see that the criterion for 
efficient reheating can be achieved for moderate values of 
the parameters. We only need for the length-scale cutoffs 
1 /A and e to be somewhat smaller than the string scale, 
while the size of the extra dimensions should exceed a 
few times / s . Using l|66l) . we can write the criterion for 
efficient reheating as 



(2na(M s )) 1/i l s = 0.7l s 



(68) 



We have not taken into account the back-reaction 
of the particle production on the tachyon background, 
which is why our calculation allows for more reheating 
than is energetically possible. The back reaction will sup- 
press somewhat the actual efficiency of reheating, but we 
don't expect a dramatic reduction. Given that we have 
been conservative in our estimates, such as ignoring the 
contributions from produced KK photons which will de- 
cay into massless photons, our result makes it plausible 
that a large fraction of the original energy can be con- 
verted into visible radiation. 



VIII. SUMMARY AND CONCLUSION 

We have argued for the possibility that our visible uni- 
verse might be a codimension-two brane left over from 
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annihilation of a D5-brane/antibrane pair at the end of 
inflation. In this picture, reheating is due to produc- 
tion of standard model particles (e.g. photons) on the 
final branes, driven by their couplings to the tachyon 
field which encodes the instability of the initial state as 
well as the vortex which represents the final brane. We 
find that reheating can be efficient, in the sense that a 
sizable fraction of the energy available from the unstable 
vacuum can be converted into visible radiation, and not 
just gravitons. 

The efficiency of reheating is greatest if the radius of 
compactification of the extra dimensions is larger than 2- 
3 times the string length l s . The efficiency also depends 
on phenomenological parameters we had to introduce by 
hand in order to cut off ultraviolet divergences in the cal- 
culated particle production rate: namely e, a nonvanish- 
ing radius for the final brane, and A, an explicit cutoff on 
the momentum of the photons produced. The latter must 
be introduced to correct for discontinous time deriva- 
tives in our simplified model of the background tachyon 
condensate; the actual behavior of the condensate corre- 
sponds to a cutoff of order A ~ X/l s . It is less obvious why 
the effective field theory treatment should give divergent 
results as the thickness of the final brane (2e) goes to zero, 
but it seems clear that a fully string-theoretic computa- 
tion would give no such divergence, and therefore it is 
reasonable to cut off the field-theory divergence at e ~ l s . 
Given only these mild assumptions, our estimates H66I67(I 
predict that the fraction of the available energy which is 
converted into visible radiation is p/p c — na(M a ) = 0.25. 
This simple estimate counts only photons; in a more re- 
alistic calculation, it would be enhanced by the number 
of light degrees of freedom which couple to the tachyon, 
which could be much greater 1. Moreover it could also 
be enhanced by the production of massive KK modes, 
which correspond to string excitations when the vortex 
has formed [SHI 

Our analysis makes significant improvements to the 
previous work of We considered the formation of a 
vortex in the tachyon field rather than a kink, in line with 
the descent relations for stable Dp-branes. We found an- 
alytic solutions for the tachyon field which give the time 
dependence in the vicinity of the defect while it is form- 
ing, for both the vortex and the kink solutions. In the 
latter case we verified that this solution reproduces the 
known dynamics of kink formation which was determined 
numerically in 10], giving us more confidence in the vor- 
tex solutions, which are quite analogous. Our explicit 
solution for the tachyon background is nevertheless too 
complicated for computing the production of particles 
on the defect. We therefore approximated it by a sim- 
pler ansatz with the same qualitative behavior, which 
allows for analytic solutions of the gauge fields in the 
background. This ansatz resembles a gauge field in a 6D 
spacetime with two compact spatial dimensions which 
are contracting with time, and leads to simple analytic 
results for the energy density of photons produced during 
the contraction, in the regime where the extra dimensions 



are large compared to the string scale. 

There are still some outstanding questions to be ad- 
dressed concerning this scenario. First, we have made 
reference to the Kibble mechanism for the creation of 
the final state defects. If we assume the causal bound 
of one defect per Hubble volume then this would imply 
that the size of the extra dimensions must exceed the 
inverse Hubble rate; otherwise there would be enough 
time for the fields to straighten themselves out and the 
putative vortex-antivortex pair would immediately anni- 
hilate. For example if we take the string scale M s to 
be 10 16 GeV then H ~ M%/M p and we would need the 
compactification scale to be of order Rq ~ (M p /M s )l s . 
Our results indicate that efficient reheating is compati- 
ble with a large compactification scale. However, taking 
the remaining four extra dimensions (which the initial 5- 
brane/antibrane pair do not wrap) to be string scale is 
not consistent with getting inflation since the initial state 
5-brane and antibrane cannot be sufficiently separated to 
satisfy the slow roll conditions. 

One the other hand, our results indicate that the re- 
heating can be efficient for Rq only a few times the string 
length, though in this scenario a naive application of the 
Kibble mechanism does not favor having the final state 
defects span the three large dimensions. These require- 
ments may not be prohibitive since the question of how 
these defects form dynamically at the end of inflation 
is a quantitative one which merits further investigation. 
In principle the correlation length for the initial fluctua- 
tions of the tachyon field could be as small as the string 
length. We point out also that it is possible that the 
dynamics of the formation of tachyon defects is qualita- 
tively different from defect formation in a conventional 
scalar field theory. For example, the numerical investi- 
gation of |10|. it was found that small kinks in the ini- 
tial configuration which are in causal contact with each 
other do not dynamically straighten themselves out as 
they would in a conventional, nontachyonic field theory. 
Instead, every place where the field crosses zero in the 
initial state develops a full-blown kink, so long as there 
was enough energy in the bulk to produce the required 
number of kinks. Another indication that the dynamics 
of the tachyon field may be qualitatively different from 
an ordinary scalar field theory comes from the [29j in 
which the causal structure of the tachyon Dirac-Born- 
Infeld action was studied. The authors of 29] found that 
small fluctuations of the tachyon field propagate accord- 
ing to an effective metric which depends on the tachyon 
background. In the case of a homogeneous rolling back- 
ground it was found that as the condensation proceeds 
the effective metric contracts to the Carroll limit of the 
Lorentz group so that the tachyon light cone collapses 
into a timelike half line and the tachyon fields at different 
spatial points are decoupled. We feel that quantitatively 
determining the dynamics of the formation of tachyon 
defects at the endpoint of D-brane inflation is a question 
which deserves further investigation. 

In the present analysis we have not included the grav- 
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brane 
stack 

FIG. 9: Warped compactification with branes localized in 
throats on opposite sides of a stack of branes. 

itational or Ramond-Ramond forces between the vortex 
and antivortex which would attract them toward each 
other and lead to their eventual annihilation. How do 
we insure that the braneworld on which we are supposed 
to live is safe from annihilation with an antibrane in the 
bulk? One possibility is to have warping caused by a 
stack of branes which wraps only the equator of the ex- 
tra dimensions, as illustrated in fig. [5] Such a braneworld 
scenario using the AdS soliton solution for the bulk has 
been considered in ref. 0, H3, EJ. The advantage for 
our scenario is that the warping can provide a barrier to 
the annihilation of the brane-antibrane pair, since it is 



energetically favorable for them to remain within their 
respective throats. 

In solving for the tachyon background we have also 
ignored the possibility of caustic formation in the bulk 
|17| by taking initial profiles without too much curva- 
ture. It is possible that caustic formation may be an 
artifact of the derivative truncation which leads to the 
Born-Infeld type of Lagranian for the tachyon. See 
for a discussion of the problems of dynamical equations 
with infinitely many derivatives. 

Brane-antibrane inflation and braneworld cosmology 
are two of the most important applications of string theo- 
retic ideas to cosmology. We find it intriguing that these 
two ideas might be combined in the way we have de- 
scribed. An outstanding challenge is to find some observ- 
able signatures that would be able to test our scenario, 
for example through the gravitational wave component 
which is expected to be a major component of the radi- 
ation produced during reheating. 
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ing the Kibble mechanism. Our research is supported by 
the Natural Sciences and Engineering Research Council 
(NSERC) of Canada and Fonds quebecois de la recherche 
sur la nature et les technologies (NATEQ). 
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